THE EXPECTED NUMBER OF ZEROS OF A 
RANDOM SYSTEM OF p-ADIC POLYNOMIALS 



STEVEN N. EVANS 

Abstract. We study the simultaneous zeros of a random family 
of d polynomials in d variables over the p-adic numbers. For a 
family of natural models, we obtain an explicit constant for the 
expected number of zeros that lie in the d-fold Cartesian product 
of the p-adic integers. Considering models in which the maximum 
degree that each variable variable appears is N, this expected value 

/Llog, (1 + p-l + + . . . + p-d) -1 

for the simplest such model. 



1. Introduction 

Various questions regarding the distribution of the number of real 
roots of a random polynomial were con sidered in [L038, L039, L043] 
and were taken up in |Kac43b[ IKac43al IKac49j , where the main result 
is that the expected number of roots of a degree n polynomial with 
independent standard Gaussian coefficients is asymptotically equiva- 
lent to I log n for large n. There has since been a huge amount of 
work on various aspects of the distribution of the roots of random 
polynomials and systems of random polynomials for a wide range of 
models with coefficients that are possibly dependent and have distri- 
butions other than Gaussian. It is impossible to survey this work 
adequately, but some of the more commonly cited early papers are 
|LS68a[ ILS68b[ llM71al lIMTlbj . Reviews of the literature can be found 



BRS86t rEK95t lEKQGt IFarQSj . and some recent papers that indicate 



m 

the level of sophistication that has been achieved in terms of results 
and methodology are pTV95l ITZ97l IBRM IbI^ IDPSZn21 ISZn3a[ I^^ZOil 

IEzo3b, wEnSj. 
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In this paper we study the roots of random polynomials over a field 
other than the real or complex numbers, the field of p-adic numbers for 
some prime p. Like the reals, the p-adics completion of the 

rationals with respect to certain metric - see below. They are the pro- 
totypical local fields (that is, non-discrete, locally compact topological 
fields) and any local field with characteristic zero is a finite algebraic 
extension of the p-adic numbers (the local fields with non-zero char- 
acteristic are finite algebraic extensions of the p-series field of Laurent 
series over the finite field with p elements). 

In order to describe our results we need to give a httle background. 
For a fuller treatment, we refer the reader to |Sch84j for an excellent 
introduction to local fields and analysis on them. 

We begin by defining the p-adic numbers. Fix a positive prime p. 
We can write any non-zero rational number r G Q\{0} uniquely as 
r = p'^{a/h) where a and h are not divisible by p. Set |r| = p~^. If we 
set |0| = 0, then the map | ■ | has the properties: 

|x| = a; = 0, 
(1) \xy\ = \x\\y\, 

\x + y\ < \x\\/ \y\. 

The map (x, ?/) ^ \x — y\ defines a metric on Q, and we denote the 
completion of Q in this metric by Qp. The field operations on Q extend 
continuously to make Qp a topological field called the p-adic numbers. 
The map | ■ | also extends continuously and the extension has properties 
(P). The closed unit ball around 0, Zp := {x E Qp : \x\ < 1}, is the 
closure in Qp of the integers Z, and is thus a ring (this is also apparent 
from ((T))), called the p-adic integers. As Zp = {x G Qp : |a;| < p}, 
the set Zp is also open. Any other ball around is of the form {x G 
Qp : |x| < p^^} = p'^Zp for some integer k. Such a ball is the closure 
of the rational numbers divisible by p^, and is thus a Zp-sub-module 
(this is again also apparent from ([Q)). In particular, such a ball is an 
additive subgroup of Qp. Arbitrary balls are translates (= cosets) of 
these closed and open subgroups. In particular, the topology of Qp has 
a base of closed and open sets, and hence Qp is totally disconnected. 
Further, each of these balls is compact, and hence Qp is also locally 
compact. 

There is a unique Borel measure A on Qp for which 

\{x + A) = \{A), xGQp, 
X{xA) = \x\X{A), X G Qp, 
A(Zp) = 1. 
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The measure A is just suitably normalized Haar measure on the additive 
group of Qp. The restriction of A to Zp is the weak limit as n — > oo 
of the sequence of probability measures that at the n-th stage assigns 
mass p"" to each of the points {0, 1, ... — 1}. 

There is a substantial literature on probability on the p-adics and 
other local fields. Two notable early papers are |M^d85t iMadQOj. We 
have shown in a sequence papers |Eva89| lEvaQH lEva93| IEva95| lEvaOTBl 
lEvaOlal IEva02j that the natural analogues on Qp of the centered Gauss- 
ian measures on R are the normalized restrictions of A to the compact 
Zp-sub-modules p'^Zp and the point mass at 0. More generally, the nat- 
ural counterparts of centered Gaussian measures for are normalized 
Haar measures on compact Zp-sub-modules. We call such probability 
measures Qp-Gaussian and say that a random variable distributed ac- 
cording to normalized Haar measure on Z^ is standard Qp-Gaussian. 
There are numerous papers on Markov processes taking values in local 
fields for exa mple fTRFH [KKM IAK/991 [KkM [KtM [KtM [KtM 
IKZ04| ISJZOSj . There are also extensive surveys of the literature in the 
books |Khr97[IKI^IKN()4j . 

If we equip the space of continuous functions C(Zp,Qp) with the 
map / 11/11 := sup{|/(t)| : t G Z^}, then || ■ || is a p-adic norm in the 
sense that 

11/11 = 0^/ = 0, 
||a/|| = |a|||/l|, aeQp, /eC(Zj,Qp), 
ll/ + ^7ll< 11/11 V||^7||. 



Moreover, C{Z'^,Qp) is a p-adic Banach space in the sense that it is 
complete with respect to the metric (/, (?) ^ \\f — g\\- 

There is a natural notion of orthogonality on the space C(Zp, Qp). A 
collection {/o,/i, . . .} is orthogonalii \\ Ylk=o^kfk\\ = VLo 
any n and any G Qp. At first glance, this looks completely unlike the 
notion of orthogonality one is familiar with in real and complex Hilbert 
spaces, but it can be seen from ^Sch84j that there are actually close 
parallels. It is apparent from |Sch84j that the sequence of functions 



{t ^ (Dlfc^O) where (*) := (the Mahler basts) is a very 

natural orthonormal basis for C(Zp,Qp) (that is, it is orthogonal and 
each element has unit norm). It is not hard to see that the functions 

{ti,t2,...,td) ^ ( !M ( M ■ ■ ■ ( M ' < h,k2, - ■ ■ ,kd < oo, 



are an orthonormal basis for C(Zp, Qp). 
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Putting all of these ingredients together, we see that a natural model 
for a random system of d independent identically distributed Qp- valued 
polynomials in d variables lying in Zp is the system 

t2, . . . , t.) := E -^Z.,, (^^J Q^J ■ ■ ■ Q^j , l<^<d, 

where the sum is over multi- indices k = {ki, k2, ■ ■ ■ , kd), for each i the 
constants G Qp are zero for all but finitely many k, and the Qp- 
valued random variables are independent and standard Qp-Gaussian 
distributed. 

Assumption 1.1. Assume that 7^ and Oe^. 0, 1 < j < d, where 
ei := (1, 0, 0, . . . , 0), 62 := (0, 1, 0, . . . , 0), and so on. By re-scaling, we 
can assume without loss of generality that Oq = 1. We will also suppose 
that |afc| > \ai\ when A; < ^ in the usual partial order on multi- indices 
(that is, if k = (fci, k2, ■ ■ ■ , k^) and £ = ^2, • • • , ^d), then kj < ij 
for 1 < i < d). It follows from the orthonormality of the products of 
Mahler basis elements that each {Fi, F2, . . . , F^) maps Z;^ into Z'^. 

Theorem 1.2. Suppose that Assumption ll . l\ holds. For {xi,X2, ■ ■ ■ , Xd) € 
Zp, the expected number of points in the set 

{(ti, t2, . . . , U) e Zj : t2, . . . , td) = x„ 1 < i < rf} 

is 

{I + p-^ + p-'^ + ■ ■ ■ + p-'^y^ . 

Following some preliminaries in Section |2l we give the proof in Sec- 
tion |21 However, we provide a heuristic argument now as motivation 
for the development we need to do in Section |21 Because we are arguing 
heuristically, we do not justify various interchanges of limits, sums and 
expectations. 

Suppose first that d = 1 and x = 0. Write F = ^fcafcZfc(^'J for 
Fi. Let Bn^o,Bn^i, . . . ,Bn,p"-i be a list of the balls of radius in 
Zp, numbered so that G -Bn.o- Let J™'" be the indicator of the event 
that the graph of F intersects Bm,i x Bnj- The number of zeros of F, 
|{t G Z : = 0}|, is 

lim lim V/™o'"- 

m— >oo n— >oo ' ' ' 

i 

Because Zq is distributed according to Haar measure on Zp, the distri- 
bution of z + F is the same for all z & Zp and so the expectation in 



nv 

17=1 h=i 



I Clhe 



h 
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question is also the expectation of 

lim lim > > /, 



m,n 

id ■ 



As we observe in Sectional -F is a stationary process on Zp (this is not 
at all obvious and will hold if and only if |ao| > |ai| > . . ., hence our 
assumption to this effect). Consequently the expectation in question is 
also the expectation of 

lim lim V I^f. 

m— >oo n— ►oo ' ^ '-^ 

j 

As in the real case, polynomials look approximately like affine func- 
tions on small scales, so for large m the restriction of F to the ball 
Bm,o of radius p"*" around is equivalent up to first order to a random 
affine function t i— > + bWt, where W is standard Qp-Gaussian and 
6 G Qp is a non-zero constant. This implies first of all that for large 
m the restriction is injective, so that lim„^oo P""" ^ j -^oj'" is the Haar 
measure of the image of -Bm,o by F. Moreover, the image has Haar 
measure approximately that of the image by t + bWt, which is 

exactly |6| ll^lp"™". Thus the expectation in question is nothing other 
than the expectation of It remains to note that \W\ takes the 

value with probability p"*" — for r = 0, 1, 2, ... to conclude 

that the expectation of \bW\ is \b\ J2ri^ — p^^)p~'^'^ = \b\{l +p~^). 

Essentially the same heuristic argument works for general d. Once 
again the problem is reduced to considering the expected Haar measure 
of the image of a small ball by a random affine function. Computing 
the actual value of the expectation is more complicated however, as it 
involves evaluating the expected value of the determinant of the linear 
part of the affine function. 

This paper appears to be the first to consider roots of random poly- 
nomials over the p-adic field. There has been some work on random 
polynomials over finite fields, see |( )do92[ IXHTOSl IIM961 IPan( )4I IDF()4j . 



2. Preliminaries 

Write Arf for the rf-fold product measure A®"'. Thus is Haar mea- 
sure on the additive group of Qp normalized so that Ad(Zp) = 1. The 
Euclidean analogue of the following result is well-known. 

Lemma 2.1. For a Borel set A C and a d x d matrix H, the set 
H{A) has Haar measure \diH{A)) = \ det{H)\\d{A) . 

Proof. If H is singular, then the range of if is a lower dimensional 
subspace of Qp and the result is obvious. 
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Suppose then that H is invertible. Write GL{d,1jp) for the space 
of d X d matrices that have entries in Zp and are invertible with the 
inverse also having entries in Zp. By Cramer's rule, a matrix W is 
in GL{d,'Lp) if and only if it has entries in Zp and |det(Vr)| = 1. 
Moreover, GL[d, Zp) is the set of linear isometries of Qp equipped with 

the metric derived from the norm \{xi, X2, ■ ■ ■ , Xd)\ = ViLi (see 
Section 3 of |Evan2j ). From the representation of H in terms of its 
elementary divisors, we have 

= f/diag(/^p'=^ . . . 

for integers ki,...,kd and matrices U,V E GL{d,Zp) (see Theorem 
3.1 of |Evan2p . Because |det(f/)| = |det(y)| = 1, it follows that 
|det(i7)| =p-(fci+-+M. 

From the uniqueness of Haar measure, o f/ and Xd ° V are both 
constant multiples of A^. Both U and V map the ball Z^ bijectively 
onto itself. Thus Xd° U = XdoV = Xd- 

Again from the uniqueness of Haar measure, Arfodiag(p'^^,p'^^, . . . ,p'^'^) 
is a constant multiple of A^. Now 

A, o diag(p'=^p'=^ . . . ,p'''){Zi) = Xd (Up^'^p^ 

d 

= JjA(/^Zp) = p-ik^+-+kd) 
= |det(/f)| = |det(//)|Ad(Z^). 

□ 

Write gl{d, Qp) for the space of d x d matrices with entries in Qp. 
We say that a function / from an open subset X of Qp into Qp is 
continuously differentiable if there exists a continuous function R : X x 
X — i> gl{d, Qp) such that f{x) — f{y) = R{x, y){x — y). This definition 
is a natural generalization of Definition 27.1 of jSch84j for the case 
d = 1. Set Jf{x) = R{x, x). 

The next result is along the lines of the Euclidean implicit function 
theorem. It follows from Lemma 12.11 and arguments similar to those 
which establish the analogous results for ci = 1 in Proposition 27.3, 
Lemma 27.4, and Theorem 27.5 of |Sch84j . 

Lemma 2.2. Suppose for some open subset X o/Q;^ that f : X ^ Q'^ 
is continuously differentiable. 

(i) // Jf{xo) is invertible for some xq G X, then, for all suffi- 
ciently small balls B containing xq, the function f restricted 
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to B is a bijection onto its image, f{B) = Jf{xo){B), and 
I det(J/(x))| = I det(J/(xo))| for x E B. In particular, 

Arf(/(5)) = |det(J/(xo))|A,(5). 

(ii) If Jf{xo) is singular for some Xq G X, then, for all sufficiently 
small balls B containing Xq, \d{f{B)) = o{\d{B)). 

The following result is an analogue of a particular instance of Fed- 
erer's co-area formula. The special case of this result for d = 1 and 
an injective function is the substitution formula in Appendix A. 7 of 

ISM. 

Proposition 2.3. Suppose for some open subset X of that f : 
X ^ Qp is continuously differentiate . Then, for any non-negative 
Borel function g ■ Qp ^ Qp, 

[ gof{x)\det{Jf{x))\\Mx)= [ g{y)i^f'\y)Udy). 

Proof. It suffices to consider the case when g is the indicator function 
of a ball C. Write 6 for the diameter of C. Put 

S := {x E X : Jf{x) is singular} 

and 

I := {x E X : Jf{x) is invertible}. 
From Lemma ESfii), Xd{f{S)) = 0, so that 

Xd{{y e Qi : f-\y) n S ^ ^}) = 

and 

/ g{y) Mf~\y) S) Xdidy) = 

= [ gof{x)\det{Jf{x))\\d{dx). 

From Lemma 12.2^ 11). we can cover the open set / with a count- 
able collection of balls Bk such that / restricted to Bk is a bijection 
onto its image, f{B) = Jf{xo){B) for some Xq G B, | det(J/(a:))| = 
|det(J/(xo))| for all x G B^, XdifiBk)) = | det(J/(xo))|A,(Sfc), and 
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diam/(_Bfc) < 6, so that g is constant on f{Bk). Hence 
g{y)W-\y)nBk)Udy) 



g{y) Mdy) 

g o /(x) I det(J/(x))| \didx) 



Summing over k gives 



9iy) n /) = J^go /(x) | det( J/(x))| Xdidx) 



and the result follows. 



□ 



3. Proof of Theorem 11.21 
For X G Zp, write N{x) for the number of points in the set 

{(ti,t2,...,td) e Z;^ : F,{tut2,...,U) =xu l<t<d}. 



Since Zi^ — {xi,X2, ■ ■ ■ , Xd) has the same distribution as Zi^, it follows 
that E[A^(-)] is constant. Also, by an extension of the argument for 
d = 1 in Theorem 9.3 of [ Eva89j (see also Theorem 8.2 of |Eva01bj ). 
the stochastic processes Fi are stationary. 
Thus, by Proposition 12. 3t 



E[N{x)] 



E 



E 



E[N{x)] Xd{dx) 
N{x) Xd{dx) 

\<lei{JF{t))\Xd{dt) 

E[\dei{JF{t))\]Xd{dt) 
E[|det(JF(0))|]. 



Now 



(0-l)(0- 2)... (0-/1 + 1) 
h\ 
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where the Wij are standard 



^p- Gaussian random variables and bj G 



is any constant with 



and so 



det(JF(0)) 



b=i 



det (Wi,)^ 



<i,j<d ■ 



From Theorem 4.1 in |Eva02j . we find, putting 



that 



E[|det(JF(0))|] 



■ d 



oo 



h=0 



TT TT 



h=0 



2h ^dJ^d+h-l 



The result then follows from a consequence of the g-binomial theorem, 
see Corollary 10.2.2 of |AAH,99j . 

Remark 3.1. (i) Suppose that a[ki,...,ka) = 1 ii ki < N for all i and is 
zero otherwise. Then \bj\ is just p*", where r = [logp A^J is the largest 
power of p that divides some integer ^ with 1 < i < N . 

(ii) Results about level sets of Euclidean processes are often obtained 
using the Kac-Rice formula. As shown in ,AW05j, result like the Kac- 
Rice formula are a consequence of Federer's co-area formula (see also 
|AT06j for an extensive discussion of this topic). It would be possible to 
derive a p-adic analogue of the Kac-Rice formula from Proposition 12.31 
and use it to prove Theorem II. 21 However, the homogeneity in "space" 
of (Fi, F2, . . . , Fd) makes this unnecessary. 

(iii) Because (Fi, F2, . . . , Fd) is stationary, its level sets are all station- 
ary point processes on with intensity the multiple of given in 
Theorem 11.21 

(iv) The requirement that the Fi are identically distributed could be 
weakened. All we actually use is that the distribution of (JF(0))ij does 
not depend on i. 



Acknowledgment: We thank two anonymous referees for suggestions 
that improved the presentation of the paper. 
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